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Thermodynamic properties such as magnetic susceptibility and specific heat have been computed 
for the Heisenberg Antiferromagnet with spatially anisotropic exchange on the kagome lattice on 
clusters up to N = 24 spins from the full spectra obtained by exact diagonalization. This approach 
is shown to provide a good represention of these thermodynamic properties above temperatures 
of about Jav/5 where J av is an average of the coupling constants. Comparison with experimental 
Volborthite data obtained by Hiroi et al [J. Phys. Soc. Jpn. 70,3377 (2001)] shows that Volborthite 
is best described by a model with nearly isotropic exchanges in spite of the significant distortion 
of the kagome lattice of magnetic sites in this compound and suggests that additional interactions 
are present. Comparison of the specific heat at low temperature raise the possibility that the 
density of states at low energy in Volborthite might be much lower than in the Heisenberg model. 
Magnetization curves under an applied field of the model are also investigated. The M = 1/3 
plateau is found to subsist in the anisotropic case and extend to lower field with increased anisotropy. 
For sufficient anisotropy, this plateau would then be observable for a field reasonably accessible to 
experiment. The absence of a plateau well below ~ 70 Teslas would further support a nearly isotropic 
model. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.40.Mg, 75.50.Ee 67.80.Jd 



I. INTRODUCTION 

Frustrated quantum magnetic insulators have received 
a large attention since many years, due to the possi- 
bility of observing unconventional behavior—. Amongst 
these, antiferromagnets with a kagome lattice of spin- 
1/2 appear to be promising candidates. However, de- 
spite numerous theoretical investigations, the nature of 
the ground-state of such spin systems remain an open 
question. In the most studied case of the spin-1/2 Heisen- 
berg model, numerical studies have concluded to the ab- 
sence of long-range magnetic (Neel) order and the pres- 
ence of an unusally large density of states at low en- 
ergjs2iMiSiSi2i2i£. But many questions such as the exis- 
tence of spontaneously broken symmetries or even the 
existence of a finite gap to magnetic excitations are not 
settled±&il. 

In recent years, two promissing experimental realiza- 
tion of a spin-1/2 kagome antiferromagnet have been 
synthczicd and studied for their magnetic properties: 
Cu 3 V 2 7 (OH) 2 -2H 2 0( VoZ&ort/wie) ^1 and reC ently 
ZnCu 3 (0H) 6 C1 2 (HerbeHsmUhite)^^^^.. Both 
compounds do not show any signs of ordering down to 
temperatures well below the exchange coupling strength 
and could be amongst the first realizations of a 2D quan- 
tum spin liquidii^i. Yet they appear to deviate some- 
what from a perfect realization of the kagome Heisenberg 
model. They differ in different ways. Herbertsmithite 
has a perfect kagome geometry but contains a signicant 
(probably intrinsic) percentage of impurity spins aris- 
ing from antisite disorde r 18 i 19 ' 20 and may also deviates 
from the Heisenberg model due to the presence of addi- 
tional interactions, possibly Dzyaloshinskii-Moriya inter- 



actions^^. In particular, the effect of the impurity spins 
which have significant interactions between themselves 
and with the spins residing on the kagome sites seems to 
be not easy to modelize. This complicates the analysis 
of the experimental data and the determination of the 
relevant model describing Herbertsmithite. By contrast, 
Volborthite has the advantage of being available with a 
very low impurity content. However, it presents a de- 
formed kagome geometry. There, the equilateral kagome 
triangles are distorted into isoceles triangles which sug- 
gest that two of the nearest neighbor exchange constants 
are different from the third. But to which extent the ex- 
change couplings differ and whether an Heisenberg model 
with such a spatial anisotropy of exchange couplings is 
sufficient to describe Volborthite remained an open ques- 
tion. 

The main purpose of this paper is to study this ques- 
tion by comparing the thermodynamic quantities mea- 
sured by Hiroi et a&2. with the results obtained from exact 
diagonalization (ED) of an anisotropic Heisenberg model 
(AHM), described by the Hamiltonian: 

Wahm = J ^2 SjSj + J' 22 SfcSj • (1) 

where the symbols note pair of nearest neighbour 
sites on the horizontal chains with exchange coupling 
J and (k,i) pair of nearest neighbour sites between the 
middle sites and sites on the chains with exchange cou- 
pling J' (see Fig. [1} . J and J' are taken antiferromag- 
netic (positive). The ratio a = J /J' will be used to 
measure the anisotropy of the couplings and the average 
Jav = ( J + 2J')/3 to set the energy scale. The classical 
AHM has a ferrimagnetic ground-state for a < 0.5 25 i 26 
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FIG. 1: (color online) Anisotropic kagome model. The 
kagome sites are show as black dots. J is the coupling con- 
stant between nearest neighbour spins on horizontal lines 
(solid lines in blue) whereas J' is the coupling constant be- 
tween a middle spin and its nearest neighbour spins (dashed 
red lines). 



which subsists for the spin-1/2 model over nearly the 
same range of values a. We shall thus limit ourselves to 
the region a > 0.5, where the ground state has zero mag- 
netization. Thermodynamic properties were computed 
from the full spectra of clusters up to N = 24 spins which 
is shown to yield results that may be considered as reli- 
able estimate of these quantities in the thermodynamic 
limit down to temperatures T ~ 0.2 J av or even below. 

We focus on the comparison of the experimental and 
numerical magnetic susceptibility x, acurately measured 
by Hiroi et aii 2 - over a large range of temperatures T, 
since it is presently the best known of the thermody- 
namic quantities. As discussed below, the comparison in- 
dicates that Volborthite is best described by the Heisen- 
berg model with at most a small anisotropy and reveals 
that additional interactions are present in Volborthite 
besides the nearest neighbour couplings. 

Besides the magnetic susceptibility, this paper reports 
the results obtained from these ED calculation for the 
magnetic heat capacity C V (T) and the magnetization 
curves M(H) under an applied magnetic field H. The 
available data for Volborthite are however still limited, 
which hinder the use of this quantities for the determi- 
nation of the model which could be relevant to describe 
Volborthite. The total heat capacity was measured by 
Hiroi et at^ but the magnetic part is somewhat uncer- 
tain, especially at high temperatures, since the lattice 
part is poorly known. Yet, as discussed below, the com- 
parison at low temperature raises the possibility that the 
density of states at low energy in Volborthite might be 
much lower than in the Heisenberg model. The magneti- 
zation of Volborthite was only measured by Hiroi et at^ 
under low fields (up to a few Tesla) which only gives ac- 
cess to the very low part of the magnetization curve. As 
shown below, the anisotropic Heisenberg model show a 



plateau at one third of the saturation of the magnetiza- 
tion which starts at a field that decrease with increasing 
anisotropy. Its location (or absence) will provide further 
information on the nature of the model which is relevant 
to Volborthite. 

This paper does not address the question of the nature 
of the ground state of the AHM which has been studied 
recently by Yavors'kii, Apcl and Evert o 25 i 26 or Wangei 
alii, using mean field approaches. The analysis of the low 
energy spectrum of the AHM and this question will be the 
subject of a forthcoming manuscript. ED results for the 
magnetic susceptibility and the heat capacity have been 
reported by Wangei aS^ but where not compared to the 
experimental data and where obtained only for a cluster 
of N = 12 spins. Wether this very small size systems 
could be sufficient for the comparison with experimental 
data was unclear. 

The paper is organized as follows. Numerical results 
obtained from the AHM model for the magnetic suscep- 
tibility and the heat capacity are compared with experi- 
ment in Sec. HTl and Sec. IIIIl repectively. The magnetiza- 
tion curves of the AHM model and the possibility of the 
detection of a magnetization plateau experimentally are 
considered in Sec. HVl Sec. fVl summarizes the results. 



II. MAGNETIC SUSCEPTIBILITY 

ED results for the magnetic susceptibility x(T) for dif- 
ferent cluster sizes and selected values of the anisotropy 
J are plotted in Fig. [5] to gether with the experimental 
data scaled with J av = 84. liV. This scaling was found by 
Hiroi et al to enable a fit of the data to the high temper- 
ature series expansion of the isotropic model for T > 2.5. 

A comparison of the x(T) obtained for J = 0.6,1,2 
at the different sizes indicate that the ED results be- 
come to be well converged to their thermodynamic val- 
ues down to temperature below T j J av ~ 0.2 once N > 18 
whereas the convergence of the TV = 12 results is slightly 
poorer. We can thus safely compare the N — 18 numeri- 
cal data for 0.6 < J < 2 with experiment over this range 
of temperature^Having assessed the convergence of the 
ED results, we examine the effect of the anisotropy of x 
and compare with experiment. 

In Fig. [5](a) one sees that the ED x(T) are unsensitive 
to the anisotropy for T/ J av > 1 -and coincide very well 
with the experimental data. However, anisotropy has a 
large effect on x(T) at lower T. Between T j J av ~ 1 and 
T/J &v ~ 0.1 — 0.2, where the ED results are subjected 
to little size effects and can be trusted as converged, one 
observes a strong increase of x with anisotropy which for 
a = 0.6 or 2 reach several times its value for a = 1 
and the experimental value^Qualitatively, this increase 
of the susceptibility can be understood as follows: (i) as 
a decreases toward 0.5, the system approachs the ferri- 
magnetic phase where x diverges at T = 0, (ii) for large 
a the middle spins become quasi disconnected from the 
chain spins and from one another, behaving as quasi free 
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FIG. 2: (color online) Magnetic susceptibility x vs temper- 
ature T scaled to J av = (J + 2J')/3 (see Eqn. flTJl). The 
curves obtained from exact diagonalization (ED) for differ- 
ent values of the anisotropy a = J/J' and number of spins 
N are shown by light coloured lines, specified by different 
symbols, as indicated in the legende. The black heavy line 
connects the experimental data of Hiroi et a£2 with T scaled 
to J av = 84.1 K. (a) includes ED results for 0.6 < a < 2 
and cluster sizes N = 12, 18, 24. (b) displays an enlarged re- 
gion of (a) for a better comparison of the ED data at small 
anisotropy with experimental data (only part of the ED data 
at large anisotropy are shown for clarity). 



spins under an applied field (until the field fully polarize 
these middle spins leading to the formation of a plateau 
at one third of the saturated magnetization-see Sec. lIV|) . 
As shown in Fig. [2] (a) such a large increase of x at low T 
does not occur in Volborthitc. This leads to conclude 
that a AHM with a large anisotropy is not appropri- 
ate to describe Volborthitc. However, as can be seen in 



FIG. 3: (color online) Specific heat C v vs temperature T. (a): 
results of exact diagonalizations (ED) -see caption of Fig. [2j- 
. (b): experimental results of Hiroi et with T scaled to 
J av = S4.1K. Cl ot is the total specific heat measured. C™ ag 
is the magnetic part derived by Hiroi et al by substracting 
from Cj ot the lattice contribution estimated from a Debye 
model. 



Fig. [21(b), a AHM with a ~ 1 is neither fully satisfactory. 
The shape of the ED x(T) curves differ somewhat from 
the experimental one. The AHM displays a bump in x(T) 
at T I J av ~ 1 absent in Volborthitc For small or moder- 
ate anisotropy the ED x(^) are slightly smaller than the 
experimental x(T) in the range 0.2 < T/J av ^ 1- In the 
AHM, the maximum of x{T) always occurs at a lower 
temperature than in Volborthitc. This reveals the pres- 
ence in Volborthitc of additional interactions between the 
spins besides those included in the AHM. 
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FIG. 4: (color online) Integrated entropy S vs temperature 
T. 



III. SPECIFIC HEAT 

ED results for the specific heat C V (T) are displayed 
in Fig. [3] (a) for the same values of the anisotropy a as 
in Fig. [21 For the sake of clarity experimental data are 
shown below in Fig. [3] (b). The ED and experimental 

results for the integrated entropy S(T) = L c v {x)/xdx 
are plotted in Fig. [3J 

The ED results for C v show little size effects and 
can be considered as (quasi)converged for T/J av > 0.2 
like those obtained for the susceptibility \- I n this 
range of temperature C v exhibit a broad peak around 
T/J &v ~ 0.6 — 0.8 which increases and shift to higher T 
with increased anisotropy but only very slightly. This 
is the usual main peak, associated to large energy ex- 
citations, rather insensitive to the detail of the interac- 
tions in the AHM. However, the integrated entropy under 
this peak - S(T = oo) - S(T = 0.2J av ) (see Fig. |4j) is 
only about 1/2 of the total entropy. Thus, as already 
found in the isotropic case^ an unusually large part of 
the entropy is located at low temperature, which reveals 
a large density of states at low energy. For all values of 
the anisotropy, such a feature is seen to remain. 

The exact location of this low temperature entropy - 
i.e. wether there exist one (or more) extra peak at low T— 
is however subject to some uncertainties. The size effects 
on C v at low T are large. This make the extrapolation 
of C v in the thermodynamic limit difficult. In the case 
of the isotropic model, the finite size C v , estimated for 
clusters up to N = 36, were found to display an extra 
peak at about T/J av ~ 0.1 or below-2, associated with 
a very large number of states, which is likely to survive 
in the thermodynamic limits (and a small peak at very 
low T < 0.02J av , associated to a few states, which may 
disappear in the in the thermodynamic limit) . The exact 
location of this extra peak in the thermodynamic limit 



is not known precisely, but its presence is necessary to 
account for the large entropy at low T, unless C v in- 
creases unsually fast with T for T — > 0, e.g. as T 13 with 
(3 < l 18 i 28 . For the same reasons a low T peak seems 
likely in the anisotropic case. The finite size results in 
Fig. [3] (a) indicate that this peak shifts to lower T with 
increasing anisotropy. 

The experimental data of Hiroi et are plotted 

in Fig. [3] (b) where we show both the total heat ca- 
pacity C*°* measured and the estimated magnetic part 
C™ as . C™ ag was derived, in absence of a non mag- 
netic isomorph, by substracting the lattice contribution 
estimated from a Debye model with a temperature of 
6d = 320 A'. As pointed out by Hiroi et al this crude 
approximation of the lattice contribution affect the ac- 
curacy of C™ aff , especially at high temperature. The re- 
sulting C" lag vanishes around T ~ 60A. By contrast, as 
seen in Fig. [2l the experimental magnetic susceptibility 
is still large at T ~ 60A. Moreover, the integrated en- 
tropy up to T = 60K is only ~ 70% of the total magnetic 
entropy (see Fig. [4]). 

The comparison between C™ a9 which is only accurate 
at low T and the ED C v which are only converged at 
high T is not easy. C™ ag is quite different from the ED 
C v of Fig. H (a). C™ ag exhibits a peak at T ~ 30A. The 
height of this peak is about twice larger than the height 
of the main peak of the ED results. If the temperature 
of the experiment are scaled to J av = 84. IK, this peak 
is located at a temperature T/J av ~ 0.3 This is about 
half of the temperature of the main peak of the ED C v . 
The differences for T j J av > 0.2 are thus quite large. But 
(jmag j s ra ^i icr uncertain in this range of temperature. 
Most likely, the peak of the real C™ ag is located at a 
temperature higher than ~ 30A. Yet one may notice 
that the integrated entropy up to to 60A (see Fig. HJ) is 
similar to the one found for the AHM up to T j J av ~ 0.7. 

At lower T, however, the lattice contribution is small 
and C™ ag more accurate. But there also, one sees quite 
distinct behavior. C™ a9 increases much less rapidly with 
T. There is no low T peak. As a result the shapes of 
the integrated entropy widely differs (see Fig. HJ). The 
integrated entropy obtained from C™ ag is much smaller 
than in the AHM which reveals that the density of states 
at low energy could be much smaller in Volborthite than 
in the AHM model. This may further suggest that some 
modification of the AHM are needed to describe Volbor- 
thite. 



IV. MAGNETIZATION CURVES 

The magnetization of the isotropic model a = 1 has 
already been the subject of many studies. The model 
displays a plateau at M = m/m s — 1/3 of the saturated 
magnetization m Sl which was interpreted as valence bond 
state with a y3 x y/3 superstructure^ and a jump of 
height SM = 2/9 to saturation which arises degenerate 
localized magnons^. The M = 1/3 plateau was aslo 
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FIG. 5: (color online) Magnetization curves for different value 
of the anisotropy a, obtained from exact diagonalization on 
clusters of N = 18 (dotted red lines) and N = 24 (dashed 
blue lines) spins. The magnetization M is normalized to 1. 
H s the saturation field. 



reported to subsist in the anisotropic case^i. But this 
case remained much less studied. 

The magnetization curves obtained from the ED spec- 
tra are ploted in Fig. [5] for different value of the 
anisotropy a in function of H/H s where H s is the sat- 
uration field. The M = 1/3 plateau subsist in the 
anisotropic case. The magnetization jump close to satu- 
ration disappear for a < 1 and is replaced for a > 1 by 
a smaller jump (proportional to the number of horizon- 
tal lines in the cluster) which will disappear in the ther- 
modynamic limit. The spin structure in the M = 1/3 
plateau corresponds, for a < 1, to the ferrimagnetic 
phase (present in zero field for a < 0.5), whereas for 
a > 1 it corresponds to full polarization of the middle 
spins. One may notice that the width of the M = 1/3 
plateau increases and the lower field Hi for which it ap- 
pears decreases with increased anisotropy. For a large 
anisotropy this field become quite low. The value of 
H s is 3J' = 3J av /(2 + a) if a < 1 and 2(J' + 2 J) = 
3J av (l + 2a)/(2 + a) if a > 1 (with H s = 3J av = 3 J if 
a = 1). Fig. El indicate that H 1 ~ 0.3H S if a = 1 but 
becomes < 0.1H S for a < 0.7 or a > 2. If J av = 84K, 
H s ~ 200 Teslas, so Hi ~ 70 Teslas if a ~ 1, but for a 
large anisotropy Hi could be ~ 20 Teslas i.e. in a range 
accessible to experiments. If a compound would be well 
described by the AHM Hamiltonian, the presence of a 
large anisotropy would be signaled by an observable mag- 
netization plateau at low field, to gether with a large peak 
in the magnetic susceptibility at low tcmpcrature(as seen 
in Sec. [n}. Additional interactions to the AHM Hamil- 
tonian may have different effects on the location of the 



M = 1/3 plateau and the susceptibility But most likely 
a magnetization plateau at low field would imply a large 
peak in the magnetic susceptibility at low temperature. 
Since the latter feature is not found in Volborthite, it is 
probable that the M = 1/3 plateau (if any) is not lo- 
cated at low field but around ~ 70 Teslas. Presently, the 
magnetization has been only measured for fields up to 7 
Teslas^ 2 - (and no plateau has been seen). Measurements 
at higher fields may be worth. The absence of a plateau 
up to ~ 70 Teslas would confirm a small anisotropy. If a 
plateau is observed below, its location would be a valu- 
able information for the determination of the model ap- 
propriate to Volborthite. 



V. SUMMARY 

In this work, we have investigated the thermodynamic 
properties of the spin- 1/2 spatially anisotropic Heisen- 
berg model on the kagome lattice by means of exact di- 
agonalization in order to compare with available experi- 
mental data for Volborthite. 

The exact diagonalization results for the magnetic sus- 
ceptibility and specific heat are found to be well con- 
verged to their values in the thermodynamic limit down 
to temperatures ~ 0.2 J av (and perhaps below) for clus- 
ters of N — 24 spins. The range of temperature is signif- 
icantly larger than the one accessible with present high 
temperature series expansion calculatio n 2 ^ 32 and also 
somewhat larger than with other methods such as the 
linked cluster expansion of Ref. 22 i 23 . 

The comparison of the computed magnetic suscep- 
tibility with the experimental data show that Volbor- 
thite is best described by an Heisenberg model with at 
most a weak anisotropy. This suggests that the spatial 
anisotropy of the couplings in Volborthite is weak in spite 
of the significant distortion of the kagome lattice in this 
compound. It also reveals the necessity to introduce in 
the present model other interactions in addition to the 
nearest neighbor exchanges to allow a good fit to the 
magnetic susceptibility of Volborthite below ~ J av . 

The comparison of the numerical specific heat with ex- 
perimental data reveal quite large differences but is un- 
easy since the latter are rather uncertain, especially at 
high temperature. At low temperature, the integrated 
entropy of Volborthite appears to be much smaller than 
in the model. This suggests that the high density of 
states at low energy, which is a carateristic of the Heisen- 
berg model may not be present in Volborthite. Although, 
more accurate experimental data will be necessary to 
check this assumption, this may indicate some signifi- 
cant differences between the physics of the Heisenberg 
model and Volborthite at low energy. In Herbertsmithitc, 
where the magnetic specific heat is also uncertain, a sim- 
ilar comparison of the integrated entropy at low tem- 
perature leads to an analogous conclusion^. For both 
compounds, this point seems worth further experimental 
investigation and will put constraint on the additional in- 
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teractions which are to be introduced in the Heiscnbcrg 
model in order to better describe these compounds. 

An other carateristic property of the kagomc Hcisen- 
berg model is the presence of a magnetization plateau 
at one third of the saturated magnetization. Because 
of the large value of the average exchange in Volborthite 
(~ 84K) , the saturation field may be larger than 200 Tes- 
las and the experimental study of the magnetization is 
uneasy For weak anisotropy the plateau may appears 
around ~ 70 Teslas and its width may not be large. 
Yet, the width of the plateau is found to increases with 
anisotropy and for large anisotropy this plateau starts at 
an applied field which may reach very small value more 
accessible to experiment. In view of the experimental 
data for the susceptibility, the detection of a plateau at 
low field seems however unlikely. Presently, the magneti- 
zation has been measured for fields up to 7 Teslas^ and 
no plateau has been seen. The absence (or the presence) 
of a plateau up to much higher fields would provide fur- 
ther inside into the model that can describe Volborthite. 

Although, it may not allow a perfect description of 
Volborthite, the anisotropic Heiscnbcrg model provide a 
starting point and a reference for the understanding of 



this compound (or compounds with a distorted kagomc 
geometry that would be identified in the future) and thus 
remain worth consideration. In order to get insight into 
the influence of the anisotropy on the nature of the gound 
state and the low energy excitations, exact diagonaliza- 
tion have been carried out for clusters up to N = 36 
spins which will be reported in an other paper—. Prelim- 
inary investigations of the modification to the Heiscnbcrg 
model that are required for a better description of Vol- 
borthite have been also started and will be be persued. 
First principle investigations might help to provide fur- 
ther inside into this question and more experimental data 
would be very usefull. 
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